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Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. 
The duration and hardness distribution of GRBs has two clusters!) now understood to reflect 
(atleast) two different progenitors2. Short-hard GRBs (SGRBs; Tyo <2 s) arise from compact 
binary mergers, while long-soft GRBs (LGRBs; 75; >2 s) have been attributed to the collapse 
of peculiar massive stars (collapsars)?. The discovery of SN 1998bw/GRB 9804254 marked 
the first association of a LGRB with a collapsar and AT 2017gfo5/GRB 170817A/GW1708179 
marked the first association of a SGRB with a binary neutron star merger, producing also 
gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), 
a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) 
localization regions of GRB 200826A; X-ray and radio emission further confirm that this is 
the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the 
afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with 
being the supernova (SN). Despite the GRB duration being short (rest-frame Too of 0.65 s), 
our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest 
LGRB found with an associated collapsar; it appears to sit on the brink between a successful 
and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars 


fail to produce ultra-relativistic jets. 


On August 26, 2020, at 04:29:52 UT, the Gamma-ray Burst Monitor (GBM) on board 
the Fermi Gamma-ray Space Telescope detected GRB 200826A with duration (T90) of 1.14 + 
0.13 seconds in the 50-300 keV energy rangelIn addition to Fermi (GBM trigger 620108997), 
GRB 2008264 was detected by four other Interplanetary Network (IPN) instruments (see Methods). 


The gamma-ray properties alone do not always yield an unambiguous classification. Some 
SGRBs show afterglow and host properties akin to LGRBs, e.g. Ref. ^ some LGRBs show 
no evidence for collapsars to deep limits akin to SGRBs, e.g. Ref. S Based solely on Tul 
GRB 200826A has a SGRB probability of 6596*17. Also taking into consideration the Epeak 


parameter of a Comptonized model fit to the single spectrum over the duration of the burst (see 


Extended Data Table 5), the probability that GRB 200826A is a SGRB increases to 74% (see 


Figure [I]and Methods). However, based on rest-frame energetics, GRB 200826A is not consistent 


with the SGRB population (right panel in|Extended Data Figure 1). 


Starting 4.2 hours after the GRB, we observed 180 sq. degrees of the Fermi-GBM localization 
with the Zwicky Transient Facility (ZTF Ji9 (see Fig.D). At 17 hours after the GRB, IPN triangulated 
the source to a smaller region. ZTF20abwysqy was the only candidate that passed our alert filtering 
scheme and was also inside the IPN region (see Methods). We discovered ZTF20abwysqy at a 
brightness of g = 20.86 + 0.04 mag (AB system). The previous upper limit (g > 21.3 mag at 
50) was 17.3 hours before the GRB, as part of the nominal all-sky survey mode. In addition to the 
spatial coincidence, ZTF20abwysqy was associated with a fading X-ray counterpart!!! and variable 
radio emissioni! confirming ZTF20abwysqy as the afterglow of GRB 200826A (see 


Data Table 1]Extended Data Table 3|and|Extended Data Table 4). 


ZTF20abwysqy was discovered in a galaxy with archival detections in Pan-STARRS 1 (PS 1 jl 
and Legacy Survey (LS}4 at a LS photometric redshift of z, = 0.71 + 0.14. The offset between 
the host galaxy’s centroid and the transient is 0.18 + 0.05”, corresponding to a physical distance 
of 2.09 + 0.6 kpc. We acquired a Gran Telescopio Canarias (GTC) spectrum of the galaxy; we see 
strong [OI] and [OIIT] lines at z = 0.748 (see [Extended Data Table 2] and [Extended Data Figure] 
[2). With both spectral and photometric spectral energy distribution (SED) fitting, we infer a stellar 
mass for the host galaxy of ~ 10°” Mz (see Extended Data Figure 2] Extended Data Figure 3]and 
[Extended Data Table 6]in Methods), which is near the distribution peak for LGRB IRI while below 


the median for SGRBs based on Refi 


We model the GRB afterglow using the standard synchrotron fireball model to constrain 
parameters related to the energy and geometry of the GRB central engine (see Methods). Electrons 
in the circumburst medium are accelerated by the shock wave and reach a power-law energy 
distribution characterized by the index p, N(E) o E". This results in a SED described by a series 
of broken power laws, e.g. Ref." (See Extended Data Figure 4]and Extended Data Figure 5). The 
associated isotropic kinetic energy of Exi; = 6.0134’ x 10°? erg lies in the top 596 of the Ex: iso 
distribution for SGRBs!, but within the 90% confidence range of the LGRB energy distribution! 
(see Methods). Our data can only loosely constrain the circumburst density. The upper end of the 


distribution is consistent with the values found for the median circumburst densities!? of LGRB S, 


while the lower end is more representative of SGRB fa 


In the collapsar scenario, a high-velocity stripped-envelope SN (SN Ic-BL) should follow the 
GRB detection! To test this scenario, we used the Gemini Multi-Object Spectrograph (GMOS-N; 
see Methods) to acquire r- and ?-band images of ZTF20abwysqy on three different epochs: ~28, 
7-45, and ~80 days after the GRB trigger (epoch 1, 2 and 3 respectively). Using epoch 3 as the 
reference, we undertook image subtraction using two different subtraction algorithms (see Methods 
for details). On epoch 1, our Gemini observations show evidence of a transient with an ;-band 
magnitude of 25.45+0.15 mag (with an i-band So limit of 25.9 mag, see [Extended Data Table 1); 
there is no source in the r-band observations up to a 5o limit of 25.4 mag (see Fig. B). On epoch 
2, we do not detect a source up to a 5o limit of 25.5 mag in the ?-band and 25.7 mag in the r-band. 
Thus, at a rest-frame time of ~16 days, the foreground extinction-corrected absolute magnitude of 
ZTF20abwysqy is M; — —18.0 mag. 


To understand the source of the ;-band excess, we use Markov chain Monte Carlo (MCMC) 
and full forward modeling of all multi-band observations excluding the GMOS detection to compare 
three scenarios: an afterglow only, an afterglow plus a KN, and an afterglow plus a SN. The 
KN model is based on a best-fit template to AT2017gf0~ scaled by a compilation of SGRB-KN 
candidates?! the collapsar model uses a SN1998bw template?2 with stretch and scale parameters 
drawn from a prior that is consistent with the historical GRB-SN sample?! We dismiss the 
afterglow-only and afterglow-plus-KN models because they predict ?-band flux at the time of the 


GMOS data point that is too faint by 1.6*55 mag — inconsistent with the observations at the ~50 


level (see |Extended Data Figure 6|and |Extended Data Table 7). 


We repeated the analysis of all three scenarios including the GMOS data point in order to do 
Bayesian model comparison. The Bayes factor between the afterglow-plus-KN and afterglow-only 
models is ~ 1, indicating that neither model is favored over the other. However, the Bayes factor 
between the afterglow-plus-SN and afterglow-only model is 1055, indicating that the afterglow-plus-SN 
model is strongly favored. Based on the compilation of Ref. 23l of GRBs with associated SNe, GRB 
200826A is the shortest LGRB found with an associated collapsar, with an observed Tyo of 1.13 s 
and rest-frame Tyo of 0.65 s (see Figure|l). The second shortest LGRB with a SN is GRB 040924, 
with a rest-frame Tyo of 1.29 l 


In the conventional fireball model of a LGRB, the rest-frame duration of the prompt emission, 
ty, is the difference between the duration of the activity of the central engine, te, and the time 
required to break out of the envelope, tP. The short duration of GRB 200826A suggests that te = 
tẹ, which might imply that its central engine is active only very briefly, or that the stellar envelope 
is unusually thick, as compared to other LGRBs, perhaps requiring a nonstandard progenitor23l 
Then it is natural to infer that there must be collapsars for which te < t; that fail to produce a 
fireball. For te [Ê ty, the jet may fail to clear a path for itself and remains cocooned within the 
exploding star, yet the cocoon itself may still produce a mildly relativistic shock breakout that 
manifests as a soft, quasi-thermal, low-luminosity GRB (IIGRB) like GRB 980425/SN 1998bw or 
GRB 060218/SN2006aj^6P7l For te < ty, there may be no prompt emission at all. 


GRB 2008264, then, may sit on the brink between a successful LGRB and a failed one. The 
sign of one continuous physical parameter, t.—t,, switches a collapsar-powered GRB discontinuously 
between two dominant emission mechanisms: internal shocks in a relativistic jet or mildly relativistic 
shock breakout at the surface of the exploding star. These two mechanisms correspond to two 
widely separated and disconnected regions in the hardness-duration-fluence phase space: the traditional 
border between SGRBs and LGRBs, and the exceptionally long and soft region occupied by 
IIGRBs and X-ray flashes (XRFs). The local rate of LGRBs (~1 Gpc?yr- t} is 2 orders of 
magnitude lower than the rate of IIGRBs (300 Gpc 3yr- ! 27 and 3 orders of magnitude lower than 
the rate of SNe Ic-BL (4500 Gpc 3yr- ! 28] It is understood that LGRBs, and to a lesser extent, 
IIGRBs, are collimated and relativistically beamed, suppressing the detection of events in which 
the jet is pointed away from our line of sight. However, even allowing for a beaming correction 
to their rates of 100 for LGRBs and 1-10 for IIGRBs =“ (but no beaming correction for SNe), 
LGRBs still occur at a rate (100 Gpc~*yr~') that is up to an order of magnitude lower than that 
of IIGRBs (300-3000 Gpc yr?) or that of SNe Ic-BL, suggesting that majority of the collapsars 
fail to produce an ultra-relativistic jet and instead drive a wide-angle or nearly isotropic and only 


mildly-relativistic cocoon. 


Based on indirect evidence from the host, afterglow, and gamma-ray properties, it has been 
argued that as many as 84% of bursts detected by the Neil Gehrels Swift Observatory? (Swift) and 
40% of the Fermi bursts that are nominally SGRBs (Tọ < 258) are actually misclassified LGRB E 
If this is correct, then one would expect many more short-duration GRBs with collapsars that are 
on the edge between success and failure of the jet. The discovery of GRB 200826A, an SGRB 


imposter, lends credence to this bold claim, and would suggest that the rates of such short-duration 
LGRBs is comparable to the rate of IIGRBs, up to a few hundred Gpc ?yr !. Thus, our discovery 
upholds the hypothesis that most collapsars fail to produce jets. 
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Figure 1: Gamma-ray properties of GRB 200826A in context. The peak energy based on a 
Comptonised fit, Epeak (keV), vs. the time-integrated, To; (s), for 2310 Fermi GBM GRBs. The 
data are fit with two log-normal distributions for the two GRB classes. The color of the data 
points indicate the probability with magenta being 100% SGRB and cyan being 100% LGRB. 
GRB 2008264 is surrounded by a red star with a SGRB probability of 74% (See Methods). 
Yellow squares show LGRBs with SN-bumps, yellow circles show LGRBs with spectroscopically 
confirmed SN2I black diamonds show SGRBs with claimed KN excess?2ll 
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Figure 2: Discovery of the afterglow of GRB 200826A. It was found at the position a = 
00^27"08.542*5, 6 = --34^01738.32'7* (J2000) with an uncertainty of 0.08.’ Contours in the 
left panel represent the Fermi GBM 90% (thick) and 50% (thin) credible regions from the official 
Fermi GBM localization (dark blue) and BALROG (light blue). The filled gray squares show the 
ZTF fields observed and the dark green dots are the positions of the ZTF Night 1 optical candidates. 
In the middle panel, 3c IPN triangulation is shown in pink; the four fields of the four XRT tiled 
observations are shown as light green circles and the XRT candidates are light green dots. The 
position of the afterglow is marked by the dark green reticle. The right-hand panels are centered at 
the position of ZTF20abwysqy. The cutouts, from top to bottom, show the ZTF discovery image, 
the ZTF stacked reference image, and a false color image showing the host galaxy from Legacy 
Survey DR8. In the cutouts, North is up and East to the left. 
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Figure 3: Panchromatic afterglow and collapsar confirmation. The available multiwavelength 
light curve data is over-plotted with the best fits from the afterglowpy modeling assuming an 
ISM-like environment. Detections are shown as circles with their respective error bars, and upper 
limits are shown as inverted triangles. The optical g-, r- and ¿- bands are shown in green, red and 
yellow, the XRT 1 keV data is shown in blue, the VLA in fuchsia, while the GMRT data is presented 
in pink. We show the K-corrected light curves of three well-studied GRB-SNe (SN1998bw, one 
of the brightest; SN2006aj and SN2010bh, two of the faintest) with solid black lines. The Gemini 
GMOS-N 7-band detection is shown at day 28.28, as a yellow circle, and is consistent with the 
collapsar population. We show a cutout of the 7-band ZOGY subtraction, revealing our 25.45 mag 


detection with a ZOGY corrected score of 4. 
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Methods 


1 Discovery 


Background. The traditional paradigm of classifying a GRB based only on its gamma-ray properties 
is debated!’ Although the two types of progenitors broadly map to the time duration of the 
signals, there is no clear boundary in the bimodal distribution. Some SGRBs show afterglow and 
host properties akin to LGRBs, e.g. Ref. 7134. some LGRBs show no evidence for collapsars to 
deep limits akin to SGRBs, e.g. Ref. ???6, 


In the mid 2000s, SGRB-related breakthroughs triggered by the Neil Gehrels Swift Observatory?! 
and HETE-II? included detections of SGRBs X-ray afterglows, identification of likely host galaxies??, 
and the first SGRB optical afterglow*”*'. Yet, in the first decade of optical follow-up of Swift 
SGRBs, only about 30 optical/NIR afterglows were detected!6 4? and even fewer had associated 
redshifts. On 2017 August 17, the joint detection of the binary neutron star (BNS) merger (GW170817) 
in GWs with LIGO and Virgo*, and in gamma-rays by Fermi GBM? and INTEGRAL/65, 
unequivocally confirmed BNS mergers as at least one of the mechanisms that can produce a SGRB. 
The merger illuminated the entire electromagnetic spectrum’? and the optical/NIR emission provided 


robust evidence of a radioactively powered KN??7?, 


The bulk of our knowledge about SGRBs comes from the well-studied Swift SGRB afterglow 
sample; however, the low Swift local rate of SGRB s€ has hampered the detection of a GW170817-like 
event. On the other hand, the Fermi GBM is arguably the most prolific engine of discovery 
for compact binary mergers, as it detects ~ 1 SGRB each week, which is four times the rate 
of Swifi$-9, and comparable to the rate of the LIGO and Virgo detectors when in observing 
mode (~ 1.5 per week during 03a). Unfortunately, there is a relative paucity of electromagnetic 
observations of Fermi GRBs and LIGO/Virgo GW events alike because it is challenging to pinpoint 


them within their positional uncertainties of tens to thousands of square degrees. 


In order to better understand the phenomenology of GRBs and compact binary mergers, 
wide field of view (FOV) optical instruments have looked for Fermi GRBs and LIGO/Virgo GW 
counterparts. In 2013, the Palomar Transient Factory (PTF)® used a 7 deg? camera on the Palomar 
48 inch Oschin telescope (P48) to discover the first optical afterglow of a GRB based solely 
on a Fermi GBM localization® and subsequently found afterglows of 7 other LGRBs?. Now, 


the Zwicky Transient Facility (ZTF19772), a 47 deg? camera mounted at the P48 telescope, has 
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enabled searches an order of magnitude faster in areal and volumetric survey speed and has been 
used to search the coarse error regions of Fermi GBM SGRBs”™ and LIGO/Virgo events”. 
Both PTF and ZTF have also discovered afterglow-like transients with optical, X-ray, and radio 


emission, but no gamma-ray counterpart’**?. 


Gamma-ray detections. Fermi (GBM trigger 620108997)5!, and the IPN instruments — AGILE’, 
INTEGRAL (SPI-ACS), Mars-Odyssey (HEND), Konus- Wind 9^, and AstroSat?* detected the burst. 
The GBM localization calculated by the ground software5*56 is RA = 4.7, Dec = 35.3 (J2000 
degrees, equivalent to J2000 00h 18m, 35d 17’) with a statistical uncertainty of 1.7 degrees?!. 
The BALROG localization is consistent with this®’. The final GBM localization of the GRB was 
available a few minutes after the trigger, and it enclosed 339 (63) deg? within the 95% (50%) 


credible region. 


The light curve shows a bright GRB featuring three distinct peaks in the 50-300 keV energy 
range, with with a signal-to-noise ratio (S/N) of ~100 and a Ts, duration of 1.14 + 0.13 s. This 
GRB is within the top 36% in terms of fluence (in the 50-300 keV energy band as measured over 
Tio0) from bursts recorded in the Fermi GBM catalog®. Spectral analysis was performed using 
the RMFIT 4.4.2 light curve and spectral analysis softward"] in the interval of Tọ to Tọ +1.152s 
encompassing the Tyo start time and duration. The detector selection for the analysis consisted 
of the brightest NaI detectors with angle < 60° (N6, N7, N8, N9, NB11) and BGO detector Bl. 
Nal and BGO detector data from ~8 keV to —900 keV and ~300 keV to ~40 MeV respectively 
was used for the analysis. The Band Model®*, a standard model for GRB analysis?, provided 
best fit with parameters consistent with the initial analysis”, where Epeak = 88.9 + 3.2 keV, o 
= —0.26 + 0.07, and 6 = —2.4+ 0.1. The peak flux in the 64 ms time (10 keV - 1 MeV), 
measured from Ty 41.152 s, is 64.3 + 2.1 ph cm~? s~!. The fluence over Tyo (10 keV to 1 MeV) is 
42.6 + 0.2 x 1077 erg cm ? (the fit merit of the spectral fit is 1.09). 


A time resolved spectral analysis was performed from Tọ to To 1.152 s. The data was binned 
by S/N, with at least 20c for the S/N in each bin. For 6 of the bins, the best-fit parameters of the 
Band Model are displayed in [Extended Data Table 5| The parameters for the Band model in the 
initial time interval were poorly constrained so a cutoff power law (Comptonised) model?! was 


fitted to this interval. The Band Model was primarily used as this was the best constrained model 


“Available at the Fermi Science Support Center: https: //fermi.gsfc.nasa.gov/ssc/data/p7rep/ 
analysis/rmfit/ 
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for the time integrated analysis. The Castor C-stat statistic?” was used for its robustness when 
dealing with low count statistics, while the background was modeled as a first order polynomial 
over all qualifying detectors, as discussed above, using TTE data products. Spectral evolution of 
the Fea; parameter is observed over the duration of the burst, and is found to exhibit hardness-intensity 


tracking behavior?)”. 


The distribution of GRB 79, durations demonstrates a bimodality that is best fit with double 
log-normal components. This suggests two distinct but overlapping progenitor probabilities, and 
by extension, classifications for short and long GRBs. Using only the 759, and the parameters 
derived for the Fermi GBM distribution of bursts in RefP! we determine a SGRB probability of 
64.94%11;39 for GRB 200826A. 


According to BATSE?,, the overlap of GRB To duration distributions occurs at ~2s. In 
the past, GRB hardness ratios, quantifiable for the majority of GRBs in the absence of location 
information, have been used as a measure for the spectral hardness. Hardness ratio is presented for 
the GBM catalogue events in Ref. ©, where the data shows an anti-correlation between hardness 
and duration; however, there is a large scatter in the data. We derive a Too value of 4.2s as 
the duration where an event has equal probability of being in the short or long class? based 
on 2353 GRBs. We choose the Epeak parameter instead of hardness for separating LGRB and 
SGRB populations as it does not require set energy bands for a comparison to be made”. In 
Figure|I]we presents a log-log plot of the Too and Epeax parameter, for which we fit two log-normal 
distributions, one for each class, and the probability that an event is an SGRB is translated into a 
colour. From this analysis, the probability that GRB 200826A belongs to the short class is 7496 ?5. 


GRB 2008264 triggered Konus- Wind (KW) at Tow = 16195.106 s UT (04:29:55.106). The 
propagation delay from Earth to Wind is 2.540 s for this GRB; correcting for this factor, the KW 
trigger time corresponds to the Earth-crossing time 16192.566 s UT (04:29:52.566). The burst 
light curve shows a multi-peaked pulse which starts at about Toxw — 0.10 s and has T1oo = 0.97 s, 
To) = 0.67*013 s, and To = 0.281895 s measured in the 80-1300 keV band. Considering only 
the duration, the 75; and To of GRB 2008264 are consistent with the short GRB population in 
the KW sample. The 759 and 75, durations at which a KW-detected GRB has an equal probability 
of being short- or long-duration are 2 s and 0.7 s, respectively °°. The spectral lag between the 
20-80 keV and 80-330 keV 16 ms light curves is 30 +11 ms, consistent with the bulk of KW short 
GRBs ”. 
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During the burst, KW measured five spectra in the 20 keV-10 MeV band. The first four 
with 64 ms accumulation time cover the interval from Toxw to Tokw + 0.256 s and the fifth, 
from Tokw + 0.256 s to Toxw + 8.448 s. The time-averaged spectrum of the burst (measured 
from Tokw S to Toxw + 8.448 s) is best fit in the 20 keV—2 MeV range by the Band function 
with a low-energy photon index a = 1.26*135, a high-energy photon index 9 = —2.32*0 12, and 
a spectrum peak energy Epeak = 6711? keV (x?/d.o.f. = 50/59). The burst had a fluence of 
4.60*0 5. x 1079 erg cm~?, and a 16-ms peak flux, measured from Tow + 0.544 s, of 9.81*1 51 x 
10-9 erg cm ? s^! (both in the 20 keV-10 MeV energy range). Using the spectroscopic redshift 
z = 0.7481? we have estimated the following rest-frame parameters: the isotropic energy release 
Eiso is 7.17 * 511 x 10°! erg, the peak luminosity Liso is 2.67185? x 105? erg s7}, and the rest-frame 


peak energy of the time-integrated spectrum, FEpeak,z 18 WA keV. 


With these values, GRB 2008264 is within the softest ~ 1% of KW short GRBs in terms 
of the observed Epeak,z and is within the 1e prediction band of both the ‘Amati’ and *Yonetoku' 
relations based on 315 long/soft (Type II) GRBs with known z??-!?! , Furthermore, in the Eiso-Epeak,z 
plane, GRB 2008264 is inconsistent with the short-hard (Type I) GRB population (see [Extended] 
[Data Figure I). Thus, despite the short duration, the KW parameters of GRB 200826A imply that 
it belongs to the long/soft GRB population. 


In addition to the IPN detections, the Cadmium Zinc Telluride Imager (CZTI!”) on board 
AstroSat also detected the burst. We reanalysed the data, combining data from all four quadrants 
to create 20-200 keV light curves with 0.05 s, 0.1 s and 0.2 s bins. All light curves show a 
single pulse, with hints of sub-structure in during the rise in the smaller time bins (see 
[Data Figure 1]. We process the 0.05 s light curves with the CIFT pipeline!™, which incorporates 
better data analysis as compared to the quick-look pipeline and produces more robust results. Our 
reanalysis yields a peak time of UT 04:29:52.95 - consistent with the Fermi peak. The new value of 
Too is 0.94*0 12 s, significantly shorter than the quick-look values reported in Ref. 55, but consistent 


with Fermi and Konus-Wind. 


Optical. We ingested the GBM localization map into the GROWTH Target-of-Opportunity (ToO) 
marshal, an interactive tool design to plan and schedule ToO observations for ZTF”. The observation 
plan generated by the ToO marshal relies on gwemopt !9^196, a code that optimizes the telescope 
scheduling process for skymaps with a healpix format, like the Fermi-GBM maps. The gwemopt 


procedure involves slicing the skymap into predefined tiles of the size and shape of the ZTF 
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field-of-view, determining which fields have the highest enclosed probability, and optimizing observations 
based on airmass and visibility windows. For this purpose, we used a modified version of the 
“greedy” algorithm described in Ref. '°’, implemented within gwemopt. The resulting optimized 

plan for the optical follow-up of GRB 200826A consisted of four primary ZTF fields and one 
secondary field?| The fields were observed in the r- and g-band for 300s each, starting 4.9 hours 

after the GBM detection. The observing plan for the first night covered 186 deg”, corresponding 

to 77% of the GBM region (see FigD). Once the Konus-Wind data became available on ground at 

about 17:55 UT, the 288 arcmin? IPN error box was derived using Konus, GBM, and HEND data, 
which allowed the X-ray Telescope (XRT!*) on board Swift to initiate a ToO at about 20:45 UT. 

The IPN box was published later at about 21:30 UT!™. 


Reference images of the fields are then subtracted from the ZTF ToO observations and any 


72,110 that contains relevant information about 


high significance difference (750) generates an alert 
the transient. We queried the stream of alerts using three different tools: the GROWTH marshal!!!, 
the Kowalski infrastructure! T] and AMPEL!!3-1144] Our filtering scheme has been described in 
previous SGRB and GW searches”*-”°, but we summarize the main points here. We aim to identify 
sources that 1) are spatially coincident with the skymap, 2) are detected only after the GBM trigger, 
3) are far from known bright sources, 4) are spatially distinct from Pan-STARRS 1 (PS 13) stars 
(based on Ref. !!8), 5) have at least two detections separated by at least 15 minutes to avoid moving 
objects, 6) have a real-bogus score (RB!!?.!?) greater than 0.15, and 7) showed an increase in 
their flux relative to the reference image. For each candidates passing these filters, we visually 
inspected the light curve, cross-matched with the PS1-DR2 catalog to check for previous activity, 
and cross-matched against the Wide-field Infrared Survey Explorer (WISE'^9) catalog to determine 
whether it was consistent with an active galactic nucleus (AGN) based on its position in the WISE 
color space?'*?, Additionally, we searched the Minor Planet Center (MPC) to ensure that our 
potential counterparts were not consistent with known solar system objects. Forced photometry ^ 
was performed using data from the 10 nights previous to the trigger, to remove young SNe from the 
sample. A final quality check was done by querying for alerts around each candidate, and rejecting 
transients that have multiple alerts within a radius of 15" as these alerts could suggest artifacts, 


ghosts, or slow-moving objects near the stationary points. 


PThe ZTF secondary fields are strategically located to cover the chip-gaps of the fields in the primary grid 


"https://github.com/dmitryduev/kowalski 


“https://github.com/AmpelProject 
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After the first night of observations 28195 alerts were generated in the region and 14 sources 
passed every stage of our filtering criteria6. The SGRB afterglow, ZTF20abwysqy!”’, was not in 
our first selection as there was a ZTF alert 11" NE; however, it passed all the other filtering criteria. 
Two of the ZTF ToO fields triggered during the first night covered > 99% of the triangulated IPN 
region. However, none of the candidates reported in Ref. !*° fell there. For our second night of 


observations, we scheduled 600 s observations in r- and g-band and queried the database. 


Probability of Chance Coincidence. We roughly calculate the probability of chance coincidence 
(p) of this optical transient to be independent of GRB 2008264. For this we follow Ref. ?? and use 


Poissonian statistics to derive the probability for one or more events to be randomly coincident: 


p-1-e^? 
3 


where A is the product of three different values, À = lI ri. For this study, we consider r; to be the 


time window between the last ZTF non-detection (0.72 days before the trigger time) and the time 
of the detection (0.21 days, see[Extended Data Table T), resulting in rı = 0.93 days. The parameter 
ry is the rate of Fermi GRBs, which from the latest GBM catalog? gives ~ 0.65 bursts per day. 
Finally rs is the ratio between the total IPN area (288 arcmin?) and the sky. This derivation gives a 
p — 1.16e — 6, and allows us to rule out a random association between the GRB and the afterglow 
at the 4.870 level. 


2 Follow-up 


The summary of the follow-up results can be found in|Extended Data Table 1 [Extended Data Table 
and [Extended Data Table 4 


Optical/NIR. Here we present our optical and near infrared (NIR) follow-up results. In addition 
to our observations, MASTER? and Kitab follow-up with a clear filter??:?! led to upper limits 


of 18.3 mag and 20.4 mag respectively. 


Las Cumbres Observatory. We performed follow-up of ZTF20abwysqy with the Spectral 
camera mounted at the 2-m Faulkes Telescope North (FTN), located at Haleakala Observatory. 
Starting on 2020-08-28 11:32:26 UT we acquired three sets of images with 300s exposures each 
in g- and r-bands through the LCO observation portal] Furthermore, we obtained a second 


^https://observe.lco.global/ 


21 


(reference) epoch in g- and r-bands on 2020-09-02 UT. Our images were reduced by an automatic 
subtraction pipeline. Our pipeline retrieves images from LCO, stacks them, extracts sources from 


the image using Source Extractor! and performs photometric calibration of sources usin 
g g p p g 


the PS1 catalog. Then, the pipeline performs image subtraction using the High Order Transform 
of Psf ANd Template Subtraction code (HOTPANTS!°°) to subtract a PSF-scaled reference image 
aligned using SCAMP'**. Our photometry measurements reported in g- and r-band were determined 


after host subtraction using the LCO reference images in the same filters acquired on 2020-09-05. 


Hale Telescope. We obtained two epochs of dithered J-band imaging with the Wide Field 
Infrared Camera (WIRC!*) on the Palomar Hale 200-in (P200) telescope on 2020-08-28 UT and 
on 2020-09-04 UT, spending an hour integrating on target during each epoch. We reduced the 
images using the image reduction pipeline described in Ref. !°°. Images were aligned and stacked 
using SWarp and calibrated against the 2MASS catalog. Image subtraction between the two 
epochs was performed using the method described in Ref. '** to derive flux measurements and its 


uncertainty for the first epoch. 


Gran Telescopio Canarias. The 10.4m Gran Telescopio Canarias (GTC) (+OSIRIS) obtained 
spectroscopic observations starting on August 30, 04:30 UT. Two 1200s exposures were gathered 
with the R1000B grism and one 1200s exposure was obtained with the R2500I grism, in order to 
cover the entire 3700-10000 À range. The slit was placed covering the position of the potential 
host galaxy. Standard routines from the Image Reduction and Analysis Facility (IRAF) were used 


to reduce the data. 


Lowell Discovery Telescope. We used the Large Monolithic Imager (LMI'?) mounted on 
the 4.3m Lowell Discovery Telescope (LDT) to observe the optical transient ZTF20abwysqy on 
three different nights: August 29, September 13, and September 19 (3.2, 18.1 and 24.2 days after 
the GRB trigger). Observations were conducted with an average airmass of 1.0 while the seeing 
varied from 1.1" for the first night to 1.4" and 1.6" for the second and third night, respectively. We 
took 8 exposures of 180 s in the r-band on August 29. Images were taken with different filters 
on September 13: 5 exposures of 180 s in u-band, 4 exposures of 180 s in g-band, 6 exposures of 
180 s in ?-band and 6 exposures of 180 s in z-band. Finally, 10 exposure of 150 s in the i-band 
and r-band were taken during the last night of observations on September 19. We used standard 
procedures to perform bias and flat-field correction. The astrometry was calibrated against the 
SDSS catalog (release DR16!^?) and frames were aligned using SCAMP and stacked with SWarp. 
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After stacking the images, we extracted sources using Source Extractor and the magnitudes 


were calibrated against 45 PS1 stars in the field on average. We used HOTPANTS to perform image 
subtraction between the first and third epochs and found a source in the r-band at the location of 
ZTF20abwysky with a magnitude of 24.46+0.12 mag (see [Extended Data Table I). We determine 
the magnitudes of the host galaxy using the second epoch of observations. To verify this result, we 
used the 80 days i-band GMOS-North (see description in the paragraph below) as a reference and 
the HOTPANTS subtraction shows a source with r-band magnitude of 24.76+0.23 mag, consistent 


with the result using the LDT reference. 


Gemini Observatory. We acquired images of the transient location on September 23, October 
10, and November 7 corresponding to 28.28, 46.15 and 80.23 days after the GRB trigger respectively, 
which we denote as epoch 1, 2, and 3. We used GMOS-North, mounted on the Gemini North 
8-meter telescope on Mauna Kea, under the approved Director's Discretionary Time (DDT) proposal 
DD-104 (P.I: L. Singer). The host galaxy coordinates, in the Gemini images, are a = 00¢27'"08.5557°, 
6 = +3440138.634°. 


The first set of observations (epoch 1) was scheduled to be closest in time to a SN1998bw-like 
peak without suffering from moon illumination. All three sets of observations consisted of 14 200s 
r- and ?-band exposures, with a position angle of 45° to avoid blooming from neighboring bright 
stars. The average airmass for the observations was ~ 1.0 and the seeing was stable throughout 
the three epochs, at 0.7”, 0.5”, and 0.85”. The images were later reduced using DRAGONS/^!, a 
Python-based data reduction platform provided by the Gemini Observatory. 


We extracted sources using Source Extractor and determined a photometric zero point 


using 23 PSI stars in the field. Using the ZOGY'*” algorithm-based python pipeline, we performed 


image differencing on Gemini data. The pipeline makes use of Source Extractor,PSFEx!?, 


SWarp, SCAMP, and PyZOGY' to perform image subtraction. We took a 1300 pix x 1300 pix 
sized cutout for both the science and reference images, centered at the ZTF20abzwysqy position in 
the image in order to achieve good subtraction quality, as the images were affected by background 
variation. Sources were extracted using Source Extractor from both science and reference 
images with a 5-c detection threshold. The resulting catalogues were fed to SCAMP, which 
calculates and corrects for the astrometric errors in both science and reference catalogues with 


the help of a Gaia data release 2 catalogue!? for the same field. Using SWarp, images are 


'"https://dragons.readthedocs.io/ 
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then re-sampled to subtract the background and the point-spread function (PSF) of the images. 
Source Extractor-generated weight maps were used as input to the SWarp to generate 
variances, which when added with the images’ Poisson noise in quadrature, which results in the rms 
image. We select good sources from Source Extractor catalogues based on their S/N, full 
width at half-maximum (FWHM) and Source Extractor flag values. Sources which raised 


a Source Extractor flag were discarded. These sources were used to match the flux level of 


the images. We used PSFEx to extract the PSF model of the re-sampled images with the Source 
Extractor catalogue as input to PSFEx. The rms images, re-sampled images, PSF models and 
astrometric uncertainty are used as input to PyZOGY, which generate the difference image and 
a corrected score image as final products!?. A source is detected with a ZOGY corrected score 
of 4 in the ?-band, at the location of the transient. With ZOGY, we derive an i-band magnitude of 
25.49+0.12 mag for epoch 1. No source is detected in the r-band of epoch 1 nor in any band in 
epoch 2. 


We have confirmed the ZOGY-based results using an independent image subtraction pipeline, 
FPipe'^6, which is based on empirically measuring the PSFs of the science and reference images 
and matching them using the common PSF method (CPM'^). We detect the source with an 
i-band brightness of 25.45+15 mag in epoch 1. The results agree within uncertainties with the 
ZOGY-based subtractions. We do not detect the source in the r-band up to a 5-c limit of r > 
25.6 mag for epoch 1. No source is detected in either filter during epoch 2 up to a 5c upper limit 
of i > 25.4 mag and r > 25.5 mag! ^? (see Figure B). 


Therefore, using two independent image subtraction pipelines, we confirm the detection of a 


source in the ;-band, at the location of the transient. 


UV/X-ray. Swift began observing the IPN localization region of GRB 2008264 0.7 days after the 
trigger. The last two observations were ~4 ks and triggered as ToO observations. The XRT 
data was reduced by the the online reduction pipeline F^ 99, As the hardness ratio remains 
constant within error bars, we assume a single absorbed power law spectrum, with a Galactic 
neutral hydrogen column?! of to 6.02 x 10?°cm~?. We convert the count rates (see|Extended Data] 
to flux density at an energy of 1 keV, using the parameters derived by the Swift pipeline: 
a photon index of Tx = 1.5*97, an intrinsic host absorption of n yint = 672? x 10??cm-?, and an 


unabsorbed counts-to-flux conversion factor of 4.27 x 1071! erg cm? ct. 


$https://www.swift.ac.uk/xrt curves/00021028/ 
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The first UV/Optical Telescope (UVOT) observations of ZTF20abwysqy was ~1.6 days after 
the burst and the detection in the white filter was associated to the underlying galaxy!!! 


Radio. For the afterglow modeling, we additionally use data reported by Ref. [12] They measure 
a flux density of ~ 40 uJy at a mean frequency of 6 GHz using the Karl G. Jansky Very Large 
Array (VLA), 2.28 days after the trigger. We observed GRB 200826A in band-5 (1050-1450 
MHz) of the upgraded Giant Metrewave Radio Telescope (uGMRT'™) on 2020-09-09 15:41:30.5 
UT and 2020-09-14 23:38:25.5 UT (^ 14.5 and ~ 19.8 days after the burst respectively) under the 
approved DDT proposal ddtC147 (P.I.: Poonam Chandra). The data was recorded with the GMRT 
Wideband Backend (GWB) correlator with a bandwidth of 400 MHz divided in 2048 channels, 
with the central frequency of 1250 MHz. The total on-source time was —75 mins with overheads 
of ~30 mins. 3C48 was used as the flux and bandpass calibrator and J0029--349 was used as the 
phase calibrator. We used standard data reduction procedures in Common Astronomy Software 
Applications (CASA)? for analysing the data. The dead antennas were first flagged by manual 
inspection and the end channels were flagged due to low gain. The automatic flagging algorithms 
incorporated into the CASA task flagdata were used to remove most of the Radio Frequency 
Interference (RFI). Any remaining corrupted data were then flagged manually. The calibrated 
data was then imaged and self-calibrated to get the final image. The synthesized beam for the final 
image was c4" x 2". We found no evidence of radio emission at the GRB position on both days. 
The 3c upper limits were 48.6 uJy/beam for the observations taken on 2020-09-09 56 and 57.4 
IJ y/beam for the observations taken on 2020-09-14. 


3 The Host Galaxy 


Our GTC spectra of the host galaxy show strong [OII] and [OIII] features at z = 0.748, in 
agreement with a previous report based on data from the Large Binocular Telescope Observatory 
(LBTO)*5, 


We modeled the GTC spectrum with the Penalized Pixel-Fitting (pPXF)' to infer the properties 
of the host galaxy's stellar populations, resulting in a fit with x2 = 1.011. The pPXF results 
confirmed the [OII] features at 3726À and [OIII] lines at 4959À and 5007À, indicators of recent 
star-formation and young, hot stars. The stellar population age derived from the weighted ages 


of the templates is 1.514232 Gyr, with evidence for a younger ~0.1 Gyr population (see Figure 


Extended Data Figure 2|and Extended Data Table 6). 
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Additionally, we corrected the magnitudes of the host galaxy (see|Extended Data ‘Table 2) 


using foreground extinction maps '°* and fed them to Prospector!’, to model the SED of the 
host galaxy. Prospector uses fsps!®!*! to generate multiple stellar populations, and fits the 
observed photometry to determine the formed mass of the galaxy, age, and intrinsic extinction 
among other parameters, using the WMAP9 cosmology!9 internally. We fitted the photometry 
to a galaxy using the Chabrier'® initial mass function (IMF), the Calzetti!®* extinction curves 
for the dust around old stars, and a star formation history (SFH) with the form of te^". Our 
results using the nested sampling dynest y!9 
of Mya = 4.64 + 1.67 x 10? Mo, a mass-weighted galactic age of tga = 1.08*525 Gyr, a 
metallicity of log(Z/Z) = —1.06*0 7, a dust extinction of Ay = 0.34705, and a star formation 
rate (SFR) of 4.013':3" Mo yr! (see the resulting SED in and the 
posterior probability distributions for the free parameters in Extended Data Figure 3). We derived 


the SFR and weighted-mass age following similar studies on GRB host galaxies!667165, and the 


algorithm, gave a galactic stellar mass distribution 


stellar mass from the mass fraction derived with Prospector. 


4 Modeling 


For the modeling of the multiwavelength emission, all optical and NIR observations were from 


difference images and were corrected for foreground extinction !°°. 


The afterglow. In the standard synchrotron fireball model a power-law energy distribution characterized 
by the index p, N(E) x E™?, results in a SED described by a series of broken power laws7!70, 
The frequencies at which the broadband SED presents its breaks are the self-absorption frequency 


Va, the synchrotron frequency Vm, and the cooling frequency ve. 


The temporal decline of both the optical and X-ray data during the first four days can be 
described by a single power law model which suggests the jet-break has not yet occurred. Our ZTF 
g-band observations are the most constraining in the optical therefore we use these to estimate the 
temporal decline rate ~, = —1.05 + 0.13. We use the 1 keV XRT data to find a, = —0.89 + 0.07. 
The similar slopes between the optical and X-ray observations suggest the location of the cooling 
break frequency, ve, lies beyond X-ray frequencies. Therefore, we estimate the spectral index at 
~1 day between the X-ray and optical as Box = —0.67 + 0.02. 


We now use a, and o; to estimate the power law index of the electron energy distribution, 


p, and to determine the circumburst density profile. In a constant density (ISM-like) medium 
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QISM = Sure) (v < ve), which gives p, = 2.44 + 0.17 and p, = 2.17 + 0.10. In the wind-like 


Scenario Qwind = cas which gives p, = 1.78 + 0.17 and p, = 1.50 + 0.10. The optical to 


X-ray spectral index Sox gives an estimate for p of 2.34 + 0.04 (6 = +), Theoretical studies 


of relativistic collisionless shocks predict p 2 2 and particularly p ~ 2.2 in the ultra-relativistic 
limit '?. Given the low values of p for the wind-like scenario we choose to assume an ISM-like 


density profile and p ~ 2.4 throughout this work. 


We note the possibility of a wind-like environment with p < 2 and a spectral break 1, 
between the optical and X-rays. Such a scenario would not change the predicted optical and 
NIR emission and it would be unable to account for the late GMOS 7-band detection without 
an additional component. However, this scenario may allow for a better fit of the afterglow X-ray 
evolution, which would decay at a rate shallower than the optical by Aa = 0.25. We do not 


consider this scenario further because of the theoretically disfavored value of p. 


Bayesian afterglow modeling. We now use Bayesian inference to analyze the X-ray, optical, 
NIR, and radio counterpart. We use two independent pipelines, one using Markov chain Monte 
Carlo (MCMC) based on the EMCEE Python package!” and one using nested sampling based on 
PyMultinest!”. The pipelines use different priors and implementations, but arrive at consistent 
results. Both utilize the afterglowpy python package!’ to estimate the physical parameters 
of the multi-wavelength afterglow. Afterglowpy is a public, open-source computational tool 
which models forward shock synchrotron emission from relativistic blast waves as a function of 
jet structure and viewing angle. Descriptions of the MCMC implementation may be found in 


Refs. 176177 and the nested sampling implementation in Ref. 29 


For this work we assume Gaussian statistics for the optical and radio data, while assuming 
Poissonian statistics (via the C-statistic!?) for the X-ray data due to low detector counts. The 
Afterglowpy model is parametrized by the isotropic kinetic energy, E; Jet collimation 
angle, 0.; viewing angle, 0,; the circumburst constant density, n; the spectral slope of the electron 
distribution, p; the fraction of energy imparted to both the electrons, €e, and to the magnetic field, 


€p, by the shock. The redshift and luminosity distance of the source are held fixed. 


Our modelling of the host indicates a small galaxy of stellar mass ~ 5 x 10? M, with 
moderate extinction Ay = 0.347025 (see $3|for details). In the MCMC implementation, we 
incorporate a Small Magellanic Cloud-like host extinction correction with total-to-selective extinction 


Ry = 2.93!? implemented with the dust-extinction software package!*'. We leave the 
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color excess E(B — V) as a free parameter, with a prior distribution computed from the Ay 
posterior found from modelling the host galaxy. This posterior is similar to the distribution of 
host extinction values observed in LGRBs with well-sampled multiband photometry'®?"'**. The 


nested sampling implementation performs no extinction correction. 


Using a top hat model for the jet structure, we perform a search over the parameter space by 
allowing all the parameters Ex iso, 0c, Ou, n, p, Ee, ep, and E(B — V) to vary with broad priors. 
We report only the MCMC results, although the nested sampling results are consistent. Given both 
the highly degenerate nature of the afterglow fitting and the low number of observations available, 
some of the parameters are not particularly constrained. One exception is the spectral slope of the 
electron distribution, for which we find p = 2.4 + 0.04, consistent with the analytical results. The 
uncertainty in the circumburst density, n = 5.5* 151? x 107? cm, includes typical ranges for both 
SGRBs and LGRBs. We find very little reddening from the host, E(B — V) = 2.5133 x 107. 


The posterior probability distributions are shown in |Extended Data Figure 5| and the parameter 
estimates are listed in|Extended Data Table 7 


Bayesian model selection. In addition to an afterglow-only model, we consider an afterglow plus 


a KN and an afterglow plus a SN. 


Kilonova. We use SEDs simulated by the multi-dimensional Monte Carlo radiative transfer 


code POSSIS!*, The simulations are performed over a grid of KN parameters: dynamical ejecta 
MS", disk wind ejecta Me, opening angle ®, and the observation angle Oops (see Ref. 
for details). We use Gaussian process regression !86!87 to interpolate the model, enabling rapid 


parameter inference. 


Supernova. The SN model starts with a bicubic spline in time and frequency, LSN (t, v), that 
interpolates a K-corrected SN1998bw templat? We apply to the template a scale factor, k, and 
a stretch factor, s, which are drawn from a bivariate normal distribution that is consistent with a 
historical sample of GRB-SNe23!!®°, The model for the observed flux density is 


1 +z)k tob 
me tobs; Vobs) = 15 ( : ’ 14 Z JVobs ] - 
v ( b b ) Andi? v (1 d z)s ( ) b 


We perform Bayesian model selection to determine which model best explains the data. We 
used the nested sampling pipeline to calculate the Bayesian evidence for each of the three models. 


The Bayes factor, or the ratio of the evidences, between the afterglow-plus-KN model and the 
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afterglow-only model is ~ 1, indicating that neither model is strongly favored over the other, 
because the KN contributes negligible flux compared to the afterglow at the time of the GMOS 
observation. The Bayes factor between the afterglow-plus-SN model and the afterglow-only model 


is ~ 1095, strongly favoring the presence of a SN. 


To better understand the source of discriminating power between models, we carry out a 
posterior predictive check. Here, we sample the posterior while excluding the GMOS 7-band data 
point and then predict its value using the rest of the data. The posterior predictive distribution 
of the AB magnitude for the GMOS 7-band detection both with and without the inclusion of a 
SN contribution is shown in The fit with the SN is consistent with the 
observation, while the fit without the SN is inconsistent at the — 5o level. Therefore, inclusion of 


the GMOS 7-band data point requires a SN component, confirming a collapsar origin. 


We compared our GMOS-N detection against extinction-corrected i-band fluxes of three 
GRB-KN candidates found in the literature!*!; using photometry from GRB 130603? and GRB 160821 193-194, 
and the compiled light curve of AT2017gfd2 "6, we correct to the redshift of GRB 200826A. We 
fit each light curve to the best-fit 2D model of GW170817 from POSSIS using SNCosmo!%” 
and extract the corresponding i-band magnitude at a rest-frame time of l6days. At z = 0.748, 


both AT2017gfo and GRB 160821 would be at Mz: —10 mag, ~8 mags fainter than our detection. 


GRB 130603B does not have enough late-time detections at the same phase for comparison. 


Our detection of a source with an extinction-corrected absolute magnitude in the ?-band with 
M; = —18.0 is consistent with the population of collapsars associated with LGRBs. A typical SN 
Ic reaches its peak magnitude 10 to 20 days post-burst, at Mg = —17.66 + 1.18 mag 8. Fig. 
shows K-corrected light curves for three well-sampled GRB-SNe: SN1998bw, one of the brightest; 
SN2006aj and SN2010bh, among the faintest. 
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e Outliers 
e GRB 


Excess Rate (s^!) 


175 180 185 190 195 200 205 210 
AstroSat Time (s) +3.36112e8 


Cumulative Rate (s^!) 


175 180 185 190 195 200 205 210 
AstroSat Time (s) +3.36112e8 


Extended Data Figure 1: The AstroSat and Konus-Wind gamma-ray detections. (/eft) Upper 

pannel panel: The de-trended light curve for GRB 200826A obtained from AstroSat CZTI data. 
We combined data from all four CZTI quadrants and binned it in 0.05 s bins. We fit and subtract 
a quadratic trend from the background to obtain zero-mean data. The shaded green region and 
corresponding green symbols denote a conservative GRB time span excluded from background 
trend estimation. Similarly red points denote outliers that are automatically flagged and rejected 
from the background estimate. Lower panel: A cumulative light curve obtained by summing the 
de-trended data, and normalised such that the median post-GRB value is 1.0. The dashed horizontal 
lines denote the 5% and 95% intensity levels. The corresponding vertical dotted black lines denote 
Tos and To5, yielding Tyo of (067 0.75 s. (right) Rest-frame energetics of 331 Konus-Wind GRBs 
(SGRB: triangles, LGRB: circles) with known redshift in the Liso—Epeakz plane, with Epeak,z the 
rest frame Epeak. The hardness-intensity (‘Amati’) relation for LGRBs is plotted with its 68% and 
90% prediction intervals (dark and light gray regions, respectively). GRB 2008264, as a red star, 
appears not to be consistent with the SGRB population. 
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Extended Data Figure 2: The Host Galaxy. (/eft) In the upper panel, we show the GTC spectrum 
of the host galaxy and the lines used to determine a redshift of 0.748. In the bottom panel, the 
photometry of the host galaxy (ugrizJ, see [Extended Data Table 2) in the AB system is presented in 
red circles. The SED model and photometry from Prospector are shown in green. (right) The 
pPXF host galaxy model results described in §3. (top) The integrated spectrum (black) overlaid 
with the best-fit spectrum (orange), which sums the contributions of stars and gas in the modeled 
galaxy. The red spectrum shows the gas contribution to the spectrum, and the blue diamonds show 
the residuals to the fit. The gas is offset by 1.59e-18 erg s^! cm ?. (bottom) The pPXF weights 


(color bar) of the different stellar population templates used to construct the best-fit galaxy. 


31 


toot Ay logiZiZo) 


log(t) 


Extended Data Figure 3: 


log(M/M o) = 9.64*031 


log(Z/Zo) = -1.06*9.57 


Av = 0.347033 


tgo = 1.084338 


log(t) = 2.322105 
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log(Z/Z o) 
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Posterior distribution of the Prospector parameters. 


The 


covariances and posterior probability distributions of the parameters for the host galaxy derived 
form the Prospector modelling described in 
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Extended Data Figure 4: SED sequence of the afterglow. The SED of our model compared to 
observations at five epochs. The cooling frequency ve is located at frequencies higher than 1 keV 
(Le. v, > 2.4 x 10!* Hz). The glitches at optical v ~ 10!6 Hz are the edge of validity of our 


dust extinction model. The SN makes a large contribution at late times. See the observations in 


Extended Data Table 1|Extended Data Table 3| and|Extended Data Table 4 


5 Data Availability 


Upon request, the corresponding author will provide data required to reproduce the figures, including 


light curves and spectra for any objects. 


6 Code Availability 


Upon request, the corresponding author will provide code (primarily in python) used to produce 


the figures. 
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Extended Data Figure 5: Posterior distribution for the a£terglowpy fit. The covariances and 
posterior probability distributions of the parameters for the top hat jet afterglow with 1998bw-like 
SN model described in The afterglow parameters are the viewing angle 0,5, (rad), on-axis 
isotropic kinetic energy Fo (erg), opening angle 0, (rad), circumburst number density no (cm ?), 
electron spectral index p, fraction of energy in accelerated electrons €e, and fraction of energy 
in magnetic field eg. We assume an SMC-like extinction curve with variable E(B — V). The 
SN model is a 1998bw template with variable stretch s and scale k. The total beaming-corrected 
kinetic energy in the jet Erot (erg), computed from the afterglow parameters, is also reported. The 
histograms denote the 14, 50, and 84 percentiles of the distributions, with blue lines marking the 


solution with maximum posterior probability density. 
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Extended Data Figure 6: Posterior predictive plot for the GMOS ;-band detection. The 
posterior of the AB magnitude estimated at the time of the ?-band data point (~ 28 days after 


the trigger) using afterglow only (blue) and afterglow-and-SN (orange) light curves are shown. 
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Julian day ôt Instrument Filter AB magnitude mag Soc Limiting magnitude A, 


2459085.8814 -1.81 P48+ZTF g - - 21.50 0.21 
2459085.9678  -1.72 P48+ZTF r - - 21.46 0.15 
2459086.9011  -0.79 P48+ZTF r - - 21.33 0.15 
2459086.9696 -0.72 P48+ZTF g - - 21.31 0.21 
2459087.8986 0.21 P48+ZTF g 20.86 0.04 22.50 0.21 
2459087.9181 0.23 P48+ZTF r 20.69 0.05 22.21 0.15 
2459087.9651 0.28 P48+ZTF g 20.95 0.16 21.23 0.21 
2459088.8340 1.15 P48+ZTF g 22.75 0.26 22.53 0.21 
2459088.9038 1.22 P48-ZTF r - - 21.30 0.15 
2459088.9761 1.29 P48+ZTF g - - 21.2 0.21 
2459089.6256 1.93 P200+WIRC J - - 21.8 0.06 
2459089.9585 2.27 FNT+LCO r - - 23.30 0.15 
2459089.9698 2.28 FNT+LCO g - - 23.41 0.21 
2459090.9200 3.23 LDT+LMI r 24.46 0.12 26.37 0.15 
2459115.9675 28.28 Gemini-GMOS i 25.45 0.15 25.9 0.11 
2459115.9675 28.28 GeminitGMOS r - - 25.4 0.15 
2459133.8039 46.11 Gemini-GMOS i - - 25.5 0.11 
2459133.8039 46.11 Gemini-GMOS r - - 25.7 0.15 


Extended Data Table 1: Afterglow panchromatic observations. Observations of the 
GRB 200826A afterglow and SN. GRB 200826A was triggered at Julian day 2459087.6874. 


Julian day ôt Instrument Filter Host AB magnitude Oma, A, 
2459089.62569 9.73 P200+WIRC J 21.11* 0.16 0.05 
2459105.80457 18.11 LMI+LDT u 23.45 0.24 0.28 
2459105.80457 18.11 LMI+LDT g 23.36 0.05 0.22 
2459105.80457 18.11 LMI+LDT r 22.86 0.18 0.15 
2459105.80457 18.11 LMI+LDT i 22.66 0.16 0.11 
2459105.80457 18.11 LMI+LDT Z 22.13 0.05 0.09 


Extended Data Table 2: Host galaxy panchromatic data. Observations of the host galaxy of 
GRB 200826A. *Magnitudes are in the Vega system. 
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Julian day ôt Instrument Energy range Count rate [s !] 


2459088.386616 — 0.70 XRT 0.3-10keV — 2.16935. x 107? 
2459088.441940 0.75 XRT 0.3-10 keV — 1.91949, x 107? 
2459088.522328 — 0.83 XRT 0.3-10 keV — 1.6691. x 107? 
2459089.437860 1.75 XRT 0.3-10keV — 5.64121, x 10-? 
2459090.285972 — 2.60 XRT 0.3-10keV — 5.36122, x 10-3 
2459094.090850 — 6.40 XRT 0.3-10 keV — 3.33077... x 10-3 
2459104.352156 16.66 XRT 0.3-10keV — 1.22995. x 10-3 


Extended Data Table 3: Afterglow X-ray detections. X-ray observations of GRB 200826A. 


Julian day ôt | Instrument Frequency Flux [erg cm ? s !] 
2459089.967407 2.28 VLA 6 GHz 40 ** 
2459102.153825 14.46 GMRT 1.256 GHz < 48.6 
2459107.485017 19.79 GMRT 1.256 GHz < 57.4 


Extended Data Table 4: Radio data. Radio observations of GRB 200826A. ** VLA data from 
Ref 


Time Bins (s) Model Amplitude Epeak a B Photon Flux Photon Fluence Energy Flux Energy Fluence Fit Merit 
[keV] [ph s^! em~?] [ph cm?] [erg s^! cm7?] [erg cm~?| 
0.000-0.180 Comp 0.340.04 114.8+6.0  -0.640.1 - 32.2+0.9 5.8+0.2 3.12-0.2x1079 0.640.02 x 1075 1.1 
0.180-0.318 Band 4.2+2.3 61.3+0.5 0.510.3 -2.70.1 41.54+1.1 11.5+0.2 3.740.2x10-® — 1.2+0.2x 1079 0.9 
0.318-0.414 Band 06+0.1 141.9410.3 -0.5+0.1 -3.00.4 61.5+1.5 17.4+0.2 7.8E0.4x10-9 — 1.9+0.3x10~° 0.9 
0.414-0.506 Band 0.70.1  145.46413.5 .-0.340.] -2.30.2 61.4+1.5 23.0+0.1 9.6t-0.4x10-6  2.8+0.4x 10-6 1.0 
0.506-0.607 Band 1.2+0.4 99.647.4 -0.0440.2 -2.5+0.2 56.441.4 28.8+0.1 6.9:-0.3x1079 — 3.54-0.3x 10-5 0.9 
0.607-0.747 Band 0.9+0.2 95.345.3 -0.240.]  -3.340.4 41.74+1.1 34.6:-0.2 4.2+0.2x10-® — 4.1--0.3x 10-9 0.9 
0.747-1.152 Band 0.64-0.5 41.24+3.9 -0.340.4 -2.70.2 15.0+0.5 40.7+40.2 0.9:-0.1x 1079 — 4.54-0.3x 1075 1.1 


Extended Data Table 5: Parameters of the Fermi-GBM fit. Fitting parameter of the Fermi-GBM 
gamma-ray spectrum of GRB 2008264. 


Line F, 

A 1077" erg s7! cm~? 
[OI] 3726 1.268 + 0.37 
[OI] 3729 2.522 + 0.37 
[OIII] 5007 8.21 35 


Extended Data Table 6: Host galaxy emission line fluxes. Fluxes derived with pPXF for the 
lines detected in the GTC spectrum of the host galaxy. 
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Parameter Unit Value 


0, [rad] 0.204053 
Foy ise [erg] D x 10°? 
0. [rad] 024 os 
n [nur] 5.54573 x 107? 
p 2.4007 
Ee 0.42403 
EB 6.441974 x 1075 
E(B — V) 2.5133 x 107? 
stretch s 0.734023 
scale k 0.281025 
EX [erg] — 1.4555 x 109! 


Extended Data Table 7: Afterglow properties. Posterior afterglowpy fit model parameters 
with an SMC extinction curve, a SN1998bw template and including the final GeminitGMOS 
detection. Uncertainties are quoted at 9096. Ex is the beamed corrected kinetic energy. See §4]for 


more details. 
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